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O ' Abstract 

Within the framework of Lorentz model for description of viscoelastic medium the 
^ ■ influence of deformational defect of the shear modulus is studied on melting of 

. ultrathin lubricant film confined between the atomically flat solid surfaces. The 

, possibility of jump-like and continuous melting is shown. Three modes of lubricant 

behavior are found, which correspond to the zero shear stress, the Hooke section of 
^ - loading diagram, and the domain of plastic flow. Transition between these modes 

^ I can take place according to mechanisms of first-order and second-order phase transi- 

O ' tions. Hysteresis of dependences of stationary stresses on strain and friction surfaces 

temperature is described. Phase kinetics of the system is investigated. It is shown 
that ratio of the relaxation times for the studied quantities influences qualitatively 
on the character of the stationary mode setting. 
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' 1 Introduction 

o ■ 

CJ 

Last years the problems of sliding friction of smooth solid surfaces with thin lubricant 
^ I film between them attract considerable attention [1]. Experimental study of atomically 

• flat surfaces of mica, separated by ultrathin lubricant film, has shown that it manifests the 

properties of solid at certain conditions [2]. Particularly, the interrupted motion (stick- 
slip) was observed inherent to dry friction [3,4]. Such boundary mode is realized, if film 
of lubricating material has less than four molecular layers. This is explained as solidifica- 
tion conditioned by the compression of walls. Subsequent jump-like melting takes place, 
when shear stress exceeds a critical value due to "shear melting effect". These films are 
characterized by a yield stress, which is a characteristic of failure in solids. 
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In the previous work [5] the approach is developed according to which the transition of 
uhrathin lubricant film from solid-like to liquid-like state takes place as a result of ther- 
modynamic and shear melting. The combined analytical description of these processes is 
carried out due to self-organization of the fields of shear stresses and strain, and also tem- 
perature of lubricant film taking into account additive noises of the indicated quantities 
[6,7] and correlated fluctuations of temperature [8,9]. The influence of additive noises on 
phase kinetics of the system is analyzed [10]. 

However, there a question was not considered about the reasons for the jump-like melt- 
ing and hysteresis, which were observed in the experiments [11]. The proposed work is 
devoted to flnding of realization conditions of these features taking into account the de- 
formational defect of the shear modulus. Our approach is based on the Lorcntz model 
for approximation of viscoelastic medium [5] at the conditions when the flrst-order phase 
transition is realized. Three stationary mode are found — two solid-like, corresponding 
to dry friction, and one liquid-like, meeting the sliding. The phase portraits of the system 
are built. The various interrupted modes of friction depending on initial conditions until 
setting of equilibrium can be realized. It is also shown that in the case, when relaxation 
time of strain much longer than corresponding times for stresses and temperature, the 
hysteresis is manifested in phase portraits. 



2 Basic equations 

In the previous work [5] on the basis of rheological description of viscoelastic medium, 
possessing heat conductivity, the system of kinetic equations is obtained, which determine 
the mutually coordinated behavior of shear stresses a and strain e, and also temperature 
T in ultrathin lubricant fllm during the process of friction between atomically flat mica 
surfaces. Let us write down these equations using the measure units: 



for variables a, T, respectively, where p is the mass density, is the speciflc heat 
capacity, Tc is the critical temperature, 770 = V — 2Tc) is the characteristic value of 
shear viscosity 77, tt = pPcy/n is the time of heat conductivity, I is the scale of heat 
conductivity, k is the coefficient of heat conductivity, is the relaxation time of strain. 

Go = 7]o/^6- 



Here the relaxation time of stress the temperature T^. of atomically flat mica surfaces 
of friction, and the function g{a) = G{a)/Go are introduced, where G{a) is the shear 
modulus of lubricant depending on the stress value: 





= -a + g{a)e, 

T,e = -£ + (T - 1)(7, 

ttT = (7; -T)-ae + a^. 



(2) 
(3) 
(4) 
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(3 — const > 0. 



(5) 



At g{a) = G/Go = const equation (2) is reduced to the Maxwell-type equation for 
description of viscoelastic medium by substituting e/T„ for de/dt. The Maxwell equation 
supposes the use of the idealized Genky model. For the dependence of stress on strain a{e) 
this model is represented by the Hooke law a — Ge at £ < and constant cr^ = Ge^ at 
£ > £m {cm, £m ^re the maximal values of shear stress and strain for the Hooke section, 
e > Em results in the viscous flow with deformation rate i = {a — am)/Tj)- Actually, 
the curve of dependence (T{e) has two regions: first one, Hookcan, has the large slope 
determined by the shear modulus G, and it is followed by the more gently sloping section 
of plastic deformation whose tilt is defined by the hardening factor < G. Obviously, 
the above picture means that the shear modulus depends on the value of stresses. For 
the account of this circumstance we use the simplest approximation (5), which describes 
the above represented transition of elastic deformation mode to the plastic one. It takes 
place at the characteristic values of shear stress Up and strain Sp. It is worth noting that at 
description of structural phase transitions of liquid- like lubricant the third-order invariants 
are present, which breaks the parity of the dependence of the synergetic potential V on the 
stress (7. Therefore in approximation (5) the linear term was used a/ap {(3 = 1), instead 
of square one {a/apY (/? = 2) [12], and dependence V{a) was not already even [5]. 

Expression (3) has the form of the corresponding Kelvin- Voight equation [5,13], which 
takes into account the dependence of the shear viscosity on the dimensionless temperature 

V = ^T^Y" (6) 

Note that jointly equations (2) and (3) represent the new rheological model, since they are 
reduced to the second-order differential equation with respect to the stress a or the strain 
£. Equation (4) is the expression for heat conductivity, which describes the transmission of 
heat from the friction surfaces to lubricant film, the effect of dissipative heating of viscous 
liquid ffowing under the action of stress, and the reversible mechanic-and-caloric effect in 
linear approximation. Equations (2) - (4) formally coincide with the synergetic Lorentz 
system [12,14], in which the shear stress acts as the order parameter, the conjugate field 
is reduced to the shear strain, and the temperature represents the control parameter. As 
is known this system is used for description of both phase thermodynamic and kinetic 
transformations. It is necessary also to note that rheological properties of lubricant film 
are investigated experimentally that allows to build a phase diagram [2]. 

Dependence (5) describes hysteresis at melting of thin layer of lubricant only in coordinates 
Te — a during realization of first-order phase transition [5]. Thus, deformation curve a{e) 
is monotonous and in the case of second-order phase transition it allows us to represent 
only continuous transition. However, experimental data testify that melting of molecularly 
thin lubricant film has jump-like character [2], though may take place according to the 
mechanism of second-order phase transition. As is shown below the description of the 
indicated feature is achieved replacing dependence g{a) in (2) by g{e) = G{e)/GQ, where 

«(^) = e + ^|(^. (7) 

At this time the value of parameter (3 > 0, determining the potential character, plays a 
key role also. 



In the work [5] the melting of ultrathin lubricant film during the process of friction between 
the atomically fiat surfaces of mica represented as a result of spontaneous appearance of 



shear stresses resulting in plastic flow. This is caused by heating of friction surfaces above 
the critical value T^o = 1 + Gq/G. The initial reason for self-organization process is the 
positive feedback of T and cr on e [see (3)], conditioned by temperature dependence of 
shear viscosity leading to its divergence. On the other hand, the negative feedback of a 
and £ on T in (4) has an important role since it provides the system stability. 

In accordance with such approach the lubricant represents a strongly viscous liquid, which 
behaves like an amorphous solid — has large effective viscosity and is still characterized 
by the yield stress [2,13]. The sohd-like state of lubricant corresponds to the shear stress 
(7 = 0, since in this case equation (2) falls out of consideration^ . Equation (3), containing 
the viscous stresses, is reduced to the Debye law representing the rapid relaxation of shear 
strain during microscopic time ^ a/c ~ 10~^^ s, where a ~ 1 nm is the lattice constant 
or the intermolecular distance and c ~ 10'^ m/s is the sound velocity. Thus, equation of 
heat conductivity (4) assumes the form of simplest expression for temperature relaxation, 
which does not contain the terms representing dissipative heating and mechanic-and- 
caloric effect of viscous hquid. 

Equation (3) describes the flow of lubricant with velocity V = Ide/dt due to action of 
appearing viscous shear stress. Moreover, in accordance with Ref. [4] in the absence of 
shear deformations the temperature mean-square displacement of molecules (atoms) is 
defined by equality (-u^) = T/Ga. The average shear displacement is found from the 
relationship (u"^) = a'^a'^/G'^. The total mean-square displacement represents the sum of 
these expressions provided that the thermal fluctuations and the stress are independent. 
Above implies that the melting of lubricant is induced both by heating and influence 
of stress generated by solid surfaces in the course of friction [4]. It is shown [15] that 
the plastic flow of lubricant layer is realized at presence of clastic stress. The action of 
shear stress results to reducing of shear modulus of lubricating material. Consequently, 
the friction force decreases with increasing velocity at the contact V = Ide/dt because 
the latter leads to the growth of the shear stress according to the Maxwell stress - strain e 
relationship: da/dt— —a/T^-\-Gde/dt. It is assumed that the film becomes more hquid-like 
and the friction force decreases with the temperature growth due to decreasing activation 
energy barrier to molecular hops. 



3 Hysteresis behavior 

Let us consider the stationary states, at which the all derivatives in equations (2) - (4) 
are equal to the zero, and values o", £, T do not change in the lubricant. Then the first of 
these equations leads to the law similar to the Hooke one: 



where parameter 9 = Q/G < \, determining the ratio of the tilts for the deformation 
curve on the plastic and the Hookean regions, coefficients ge — Q /Gq < 1, and a — Sp/eg 
are introduced. 




(8) 



^ It will be shown further that non-zero stresses interval can meet the solid-like state of lubricant 
also. 



Dependence (8) at fixed a, gg, and 6 is depictured in Fig. 1. It is apparent from here that 
two situations can be realized: at small /? the curve a{e) increases monotonically (insert 
in the figure), and at 



1 + Ve 
i-Ve 



(9) 



it becomes non-monotonous. In the first case in accordance with the figure the continuous 
melting of lubricant takes place, in the second one — lubricant melts abruptly at the 
increase of stresses to the point A, and here the transition to point B occurs. At the 
further growth of the stresses the strain increases monotonically, and lubricant is liquid- 
like. If stress is decreased now to the point C lubricant conserves liquid-like structure, and 
then it abruptly becomes solid during transition at point D. At the subsequent decrease of 
stresses lubricant is solid-like. Similar transitions are represented as the first-order phase 
transitions [16], but between the states, which are not the pure thermodynamic phases. 
For understanding of these transformations the concept of shear melting is used [4]. Note 
that such hysteresis behavior was observed in the experiments [11,15]. 

Using (8), it is possible to find the abscissas of transition points A and C: 



where the sign '-' corresponds to the point A, and sign — to the point C. From Eq. (10) 
it is seen that the length of jump at melting increases with growth a, and the difference 
Ea — sc ^ with the increase of /3. Thus, at large /3 (small a) melting and solidification 
occur at the same value of strain [ea ~ £c) practically, but at different values of stresses 
a. As well as in Refs. [5]-[10], we will accept as the order parameter the shear stress a: 
at (J > (7a lubricant is liquid-like, and ii a < ac it is solid-like. In the intermediate region 
ac < cr < aA the state of lubricant is unstable, since it may exist in both phases. 

The dependence of the stationary shear stresses (Tq, a™- on the temperature of friction 
surfaces Tg is presented in Fig. 2 at parameters of Fig. 1 (Fig. 2b meets the insert in 
Fig. 1). Apparently that it is non-monotonous, and in the interval T° < < Tco the 
two-valued section is realized inherent to first-order phase transitions. The dashed curve 
corresponds to the unstable stationary values of stresses a'", the solid curve - to the stable 
(Jq. It is worth noting that a'^{Te) meets the Hookean section of dependence (7o(£o)- 

For further consideration it is necessary to write down expression for synergetic potential 
V{a). Within the framework of adiabatic approximation t^,tt <^ t„ [5] it is possible to 
set T^E 0, ttT 0, and equations (3), (4) give 



^A,C — 



6(/3 - 1) - 2 T &V(1 - - ' b = e-'-l, 



(10) 



(J -{2- Te) 



a 



(11) 
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Te + (2 - Te) 



l + (72' 



(12) 



After substitution of expression (11) in (2) we obtain the Landau-Khalatnikov equation: 



where the synergetic potential is defined by equahty 
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Note that in the region T^a < Tg < T^c this potential does not give correct result, since 
here the model has the unstable solution, describing the decrease of stresses with growth 
of strain, and this does not take into account hysteresis. For finding the form of potential 
in the indicated region we replace (2) by equation: 



where 5 is the value of the conserved stresses. In the corresponding Landau-Khalatnikov 
equation (13) the potential is fixed by expression: 



It is seen that V depends only on cr, i.e. in the situation, where stress is conserved at the 
change of temperature, the value V remains constant. 

The dependence of the potential (14) on the stress value at the fixed temperatures of 
friction surfaces is presented in Fig. 3 corresponding to the parameters of Fig. 2a. The 
hysteresis, shown in Fig. 1, is realized in this situation. We will consider this case in more 
detail. 

Below than the critical value T° stresses are absent in lubricant (ctq = 0). The form 
of potential, shown by curve 1 in Fig. 3, corresponds to this interval of temperatures. 
Here one zero minimum is realized at (Tq = 0, thus lubricant is solid-like. At point T^. — 
a plateau appears on the dependence V{o-) (curve 2). With the further increase of 
temperature in a region < Tg < Tco the potential has the form shown by curve 3. 
Here the potential barrier appears, which separates the zero and non-zero minimums of 
potential. In connection with it the system can not come in the stable state Uq ^ 0, and 
the zero value of shear stress a is realized. The dashed curve in Fig. 2a corresponds to the 
maximum of potential, the solid curve — to its non-zero minimum. As is apparent from 
figure at Tg = Tco the jump-like increase of value a takes place, and the system turns to the 
section of dependence ao(Te) (point A!). This transition is caused by that at Te = Tco the 
maximum V{(t) disappears, and with subsequent growth of the one non-zero minimum 
of potential is realized (curve 4 in Fig. 3). The expression for the critical temperature 
is obtained from the condition dV/da — 0, where V is the synergetic potential (14): 



However, section AC is unstable, because here the stress is decreased at growth of strain 
(A, fi, C and D meet the similar points in Fig. 1). In this connection further the system 
goes on the way A' - B' (this transition is already described by potential V (16), since 
stress is conserved), and passes to the plastic flow section, which corresponds to the 
hquid-hke structure of lubricant. With the subsequent increase of temperature T^, the 
stationary value of stresses (Jq grows, and lubricant becomes less viscous, here V{a) has 
the form shown by curve 4. If now the temperature of friction surfaces is decreased, the 
lubricant is liquid- like to the value Tcc, further the stress is conserved to the critical 



T^CT = -cr (5, 



(15) 



2 



(16) 



Teo = 1 + e/ge = 1 + Go/G. 



(17) 



Te = TcA (corresponding to the point A)^ . As is seen from the insert in the picture at 
Te = TcA the system passes to the point A, because here the dependence o'o(Te) becomes 
stable. At this transition system of somehow "jump over" the potential barrier (dashed 
curve of dependence). It is related to that before the jump the system is described by 
other potential V, which is characterized by absence of barrier. Here the sohdification 
of lubricant occurs, because transition takes place on the stable part of Hookc section 
(according to Fig. 1 temperature Tg < Tca corresponds to the Hooke domain). Now, the 
system is described by potential V{a) (14) with the barrier shown by curve 3 in Fig. 3. 
With the further decrease of the temperature of friction surfaces at the point Tg — 
stresses decrease abruptly to the zero, since the barrier disappears and there is one zero 
minimum of V{a). This situation corresponds to the solid-like state of lubricant also, 
but with the zero value of stress. It is supposed that the solid-like states at (Tq = and 
ctq 7^ differ by structure, because of transition between them takes place according to the 
mechanism of the first-order phase transition. Thus, the sohd-hke structure of lubricant 
at a temperature below T° is similar to the solid state, and the solid-hke structure above 
indicated temperature has the signs of the liquid state, but lubricant behaves as a solid-like 
on the whole. The different solid-like states correspond to the different modes of friction 
in the studied system [17]. At further transition to the liquid mode of friction the viscosity 
of lubricant (6) decreases, and it flows. 

At described transitions the stationary values of stresses are conserved in intervals T^a' < 
Te < TcB' and Tca < Tg < Tcc- Obviously, that equality ctq = const is satisfied with the 
increase of temperature at TcA' < Tg < TcB', since it is necessary to give some energy to 
lubricant for melting. Conservation of stresses in the region Tca < T^ < T^c at decrease 
of the sheared surfaces temperature takes place, because for transition of lubricant in 
the solid-like state it has to return energy. It is worth noting that at the conservation 
of stresses between solid-like and liquid-like phases of lubricant, it is in the intermediate 
state with differing from them structure. 

The situation, which is given by the insert to the figure 1, is shown in Fig. 2b when 
hysteresis is not realized in coordinates ao(eo)- This corresponds to the more simple case 
described in work [5]. Below the value Tco we have the solid-like state of lubricant (ao = 0), 
at Te = Tco it melts at transition to the plastic fiow region, and with the subsequent 
increase of Tg it becomes more liquid. Then, with the decrease of temperature to the 
value T° lubricant is liquid-hke (ctq 7^ 0), at Tg = T° it sohdifies abruptly (ctq = 0). 

The distinctive feature of such behavior is that whole Hooke domain is unstable, since it 
meets the maximum of potential. There is only one type of solid-like state of lubricant 
with the zero value of stress. 

Depending on the parameters of the system there the second-order phase transition can 
be observed, when the temperature T^q is moved on the left-hand side from T°. Indeed, 
the unstable stationary values of stresses cr"* are not realized, and hysteresis disappears, 
which is characterized by a presence of T°. In this case the potential, shown by curve 
1 in Fig. 3, corresponds to the solid-like structure of lubricant, which is transformed to 
the form represented by curve 4 with the increase of temperature. Thus, the continuous 
transformation takes place, since the potential barrier is absent. Here it is also possible 

^ In general case the system has to move along the curve CA, but it is unstable. Therefore it is 
necessary to introduce the described hysteresis, and to take into account that at this transition 
the system is described by potential V'{(j) (16). 



to select two situations - when hysteresis of (7o(£o) dependence is observed, and without 
it [18]. 



4 Phase kinetics 



In accordance with the experimental data for organic lubricating materials [2,6] the relax- 
ation time of stresses at normal pressure is r^- ~ 10~^° s. Since ultrathin lubricant film has 
less than four molecular layers, temperature relaxes to the value Tg during time satisfying 
the condition: 



tt < r^, T,. (18) 

According to this we will assume in equation (4) ttT 0. Also for convenience the time 
is measured in units of Tq.. As a result, we will obtain the two-parameter system in the 
form: 



a = -a + g{e)e, 



(19) 
(20) 



where r = T^jT^. The corresponding phase portraits are presented in Figs. 4 - 6 at the 
parameters of figure 2a for different Tg. 

The phase portraits describing the behavior of lubricant in the solid-like state (Tg < T°) 
are shown in Fig. 4 for different ratios of relaxation times. 

Particularly, Fig. 4a corresponds to the value r = 0.01. The isochnes obtained at equating 
to zero of derivatives in the equations (19) and (20) arc shown by dashed curves 1 and 2, 
respectively. Thus, curve 1 meets the system parameters at which stresses do not change, 
and line 2 corresponds to the case of strain conservation. These lines intersect at the origin 
of coordinates forming the unique stationary point which is a node. It is apparent that 
phase trajectories converge to the node £), i.e. stresses relax to the zero value. Thus, at 
motion on a phase plane at the arbitrary initial conditions there are two stages: on the 
first one, the instantaneous relaxation of the system takes place to the line near to isocline 
2, on the second one — the slow motion on the indicated curve. At the first stage the 
stresses are conserved that reminds the transition between the friction modes described 
in section 2. Apparently, on the final stage of system relaxation to the stationary point 
D the values of stress and strain go out in the negative domain. This can be interpreted 
as reversible motion at which the top wall of friction moves in reverse direction. Thus, 
stresses change the sign in opposite one, because of the direction of motion changes, and 
strain becomes negative also. On the other hand, it is possible to neglect the negative 
region, considering it as has no physical meaning. We can assume that at achievement 
by strain the zero value the system abruptly goes to the origin of coordinates, and the 
equilibrium sets in [10]. 

The phase portrait shown in Fig. 4b is built for the case when relaxation times of stresses 
and strain coincide (r = 1). It is also characterized by the singular point £), representing 
the node around which the weakly pronounced oscillations of short durations are realized 
now till setting of equilibrium. Here the cases are possible when stress a at first increases, 
and then decreases, and reverse situation. It means that to that moment, when the sys- 



tem will come to the equilibrium (the origin of the coordinates), the stick-slip motion is 
possible. For example, in accordance with phase trajectories, which begin at o" = 0. the 
lubricant is solid- like at first (stresses are equal to zero), then it begins to melt (stresses 
increase), and then it solidifies again (at setting of equilibrium). In order to avoid mis- 
understanding let us turn attention to the following circumstance. According to phase 
trajectories the cases are possible when with the growth of stresses the strain decreases. 
Above for the stationary values of stresses it was interpreted as instability resulting in 
hysteresis. Now, the rapid motion is realized. It allows us to assume that system is in the 
unstable state, and oscillation mode is reahzed in which the decrease of e is possible with 
the increase of a. 

Figure 4c meets the case r = 100. Here, as well as in Fig. 4a, the two stages are selected: 
the rapid relaxation to the line near to isocline 1, and further the slow motion along it. At 
the first stage the strain changes weakly, and stresses decrease very rapidly, if their initial 
values are on the right-hand side from isocline 1, or they increase at initial a on the left- 
hand side from it. At the second stage in the top part of phase portrait the configuration 
point moves along the plastic section, and in bottom one - along the Hooke domain. 
During the way of this point between the maximum and minimum of isocline 1 stresses 
increase with decrease of strain. The motion is slow here, and therefore the hysteresis has 
to be observed similar to described in Fig. 1. This feature is not shown in phase portrait, 
and its study requires the additional analysis, because there a lot of limiting cases can 
arise, which are not considered within the framework of the offered work. 

The phase portraits in Fig. 5 are presented for the same parameters and ratios of relaxation 
times as in Fig. 4. But they are built for the temperature corresponding to the section 
of Fig. 2a, for which stable and unstable stationary values of shear stresses are realized 
(r° < Te < Tco). In this case the system potential has the form shown by curve 3 in 
Fig. 3. As well as above, fines 1 and 2 are isoclines. The phase portraits are characterized 
by three singular points: by the node D at the origin of coordinates, which describes 
the dry friction; by the saddle N corresponding to the maximum of dependence V{a) 
(the unstable stationary point); by the node O meeting the non-zero stationary stress, 
which corresponds to the unstable section AC of dependence cro(Te) shown in Fig. 2a. 
These points are given by intersections of isoclines. Depending on the initial conditions 
the system may come as a result of relaxation both to the mode of stable dry friction 
(node D) and to the above described unstable section (node O). 

At r = 0.01 there is the picture shown in Fig. 5a. Here, as well as in Fig. 4a, phase 
trajectories rapidly converge to the line near to isocline 2 from any point of phase plane at 
conservation of stresses. Further the system relaxes to the nodes D, or O, and determined 
by these points the stationary modes of friction set in. The inclination of curve, along 
which motion occurs during the second stage, depends on initial conditions. So, the system 
relaxes to the point O along the plastic section of isocline 2, to the point D - along its 
Hookean one. Note that in the course of time lubricant becomes more liquid, if < en < 
uoi and vice versa more viscous at (jj > (7o, where (jj, Uiv, and uq are the initial and the 
stationary values of stresses. In these cases the system comes to the singular point O. At 
Ui < (7n the lubricant solidifies during time, and dry friction is realized (point D). 

As is obvious from Fig. 5b, for r = 1 the different types of stick-slip motion are possible 
at setting of equilibrium values of stress and strain. 



The phase portrait is shown in Fig. 5c for r = 100, where, as well as above, there are two 
stages. Since the stationary point O is on the unstable section, for the parameters of this 
figure the hysteresis is characteristic. 

The phase portraits for the temperature domain of figure 2a, corresponding to the stable 
sliding friction (Tg > T^), are shown in Fig. 6. Here in the course of evolution the non- 
zero stationary value of shear stresses is set co 7^ meeting the minimum of synergetic 
potential V{a) (the maximum of distribution function of stresses -P(cr) over their value). 
The phase portraits are characterized by two singular points — by the saddle D at the 
origin of coordinates and the node O at the non-zero values of stress and strain, which 
are given by intersections of isoclines 1 and 2. 

At T = 0.01 the situation is observed shown in Fig. 6a. Here, as well as above, phase 
trajectories converge rapidly to the line near to isocline 2 from any point of phase plane 
at conservation of stresses. Then, the system relaxes to the non-zero value cto 7^ 0, as 
a result, the stationary sliding friction sets in. However, the curve, along which motion 
occurs during the second stage, corresponds to the plastic section of dependence cr(£:), i.e. 
system is liquid-like always, except only the cases when initial value of stresses is near 
zero (the melting occurs). Note that in the course of time the lubricant becomes more 
liquid, if ui <ooi and vice versa — more viscous at ai> uq. 

As is apparent from Fig. 6b, for r = 1 at setting of stationary sliding friction the different 
types of stick-slip motion are possible. However, at any initial conditions the system comes 
to the steady shding friction. 

The phase portrait is shown in Fig. 6c for r = 100, where, as well as in Figs. 4c and 5c the 
two above described stages are seen. The stationary point O is on the plastic flow region. 
The hysteresis is characteristic also for the parameters of this figure. 

The difference of described in Figs. 4, 5, and 6 situations is achieved due to the change 
of isochne 2 form with the variation of temperature of friction surfaces, while the 
form of isochne 1 does not depend on it. The value of /3 influences substantially on the 
phase portraits, however, its variations do not change qualitatively the character of the 
system behavior. At r ^ 1 it is necessary to introduce hysteresis [18,19], because the 
slow motion of the system occurs near the isocline 1, and the condition of stationarity 
is valid in every time moment. The hysteresis is similar to described in Fig. 1, since the 
dependence, presented in this figure, represents isocline 1 shown in phase portraits. 



5 Conclusion 

The above consideration shows that the hysteresis, realized at melting of thin lubricant 
film according to the mechanism of first-order phase transition, can be described taking 
into account the deformational defect of the shear modulus. The basic feature of such 
behavior is that lubricant melts and solidifies at the different values of shear stress, which 
acts as the order parameter. The two solid-hke states of lubricant and one liquid-hke phase 
have been found, the transition between which takes place in accordance with indicated 
hysteresis. The phase kinetics is studied, and depending on initial conditions the different 
types of stick-slip motion are predicted. The elastic and thermal parameters of lubricant 
are defined, at which as a result of this motion the sliding, or the dry friction sets in. 
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Figure captions to the work of A.V. Khomenko, LA. Lyashcnko "Hysteresis 
phenomena at ultrathin lubricant film melting in the case of first-order phase transition" 



Fig. 1. The dependence of the stationary values of the shear stresses ctq on the strain Sq 
(8) at ^ = 0.2, ge = 0.6, a = 0.7, p = 5.0 (in the insert f3 = 2.0). 

Fig. 2. The dependence of the stationary values of the shear stresses aQ, (j™ on the tem- 
perature of friction surfaces Te at parameters of Fig. 1 (Fig. 2b corresponds to the 
parameters of insert in Fig. 1). 

Fig. 3. The dependence of the synergetic potential V (14) on the shear stresses a at 
parameters of Fig. 2a, curves 1-4 correspond to the values of temperature Tg = 
0.1, 0.3125, 1.0, 2.0, respectively. 

Fig. 4. The phase portraits at parameters of Fig. 2a and Tg = 0.2: (a) tt <^ = O.OItq.; 

(b) TT<^Te = Ta; (c) Tt < = IOOt^. 

Fig. 5. The phase portraits at parameters of Fig. 2a and Te — 1.0: (a) tt <^ — O.OIto-; 

(b) tt<^t,^t^; (c) tt < = lOOr^. 

Fig. 6. The phase portraits at parameters of Fig. 2a and — 3.0: (a) tt <^ — O.OItq.; 
(b) Tr < r, = r^; (c) tt < r, = lOOr^. 
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